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SEE PROJECT FOR TESTING GRAVITY IN SPACE: 
CURRENT STATUS AND NEW ESTIMATES 
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We describe some new estimates concerning the recently proposed SEE (Satellite Energy Exchange) experiment 
for measuring the gravitational interaction parameters in space. The experiment entails precision tracking of the 
relative motion of two test bodies (a heavy "Shepherd", and a light "Particle") on board a drag- free capsule. The 
new estimates include (i) the sensitivity of Particle trajectories and G measurement to the Shepherd quadrupole 
moment uncertainties; (ii) the measurement errors of G and the strength of a putative Yukawa-type force whose 
range parameter A may be either of the order of a few metres or close to the Earth radius; (iii) a possible effect of 
the Van Allen radiation belts on the SEE experiment due to test body electric charging. 



1. Introduction 

The SEE (Satellite Energy Exchange) concept of a 
space-based gravitational experiment was suggested in 
the early 90s Q and was aimed at precisely measuring 
the gravitational interaction parameters: the gravita- 
tional constant G , possible violations of the equivalence 
principle measured by the Eotvos parameter ry, time 
variations of G , and hypothetical non-Newtonian grav- 
itational forces (parametrized by the Yukawa strength 
a and range A). Such tests are intended to fill gaps 
left by current methods of ground-based experimenta- 
tion and observation of astronomical phenomena. The 
significance of new measurements is quite evident since 
nearly all modified theories of gravity and unified the- 
ories predict some violations of the Equivalence Princi- 
ple (EP), either by deviations from the Newtonian law 
(inverse-square-law, ISL) or by composition-dependent 
(CD) gravity accelerations, due to the appearance of 
new possible massive particles (partners); time varia- 
tions of G are also generally predicted [|[ [|. 

The idea of the SEE method is to study the relative 
motion of two bodies on board a drag-free Earth satel- 
lite using horseshoe-type trajectories, previously known 
in planetary satellite astronomy: if the lighter body 
(the " Particle" ) is moving along a lower orbit than the 
heavier one (the " Shepherd" ) and approaching from be- 
hind, then the Particle almost overtakes the Shepherd, 
but it gains energy due to their gravitational interac- 
tion, passes therefore to a higher orbit and begins to lag 
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behind. The interaction phase can be studied within a 
drag-free capsule (a cylinder up to 20 m long, about 1 
m in diameter) where the Particle can loiter as long as 
10 5 seconds. It was claimed that the SEE method ex- 
ceeded in accuracy all other suggestions, at least with 
respect to G and a for A of the order of metres. Some 
design features were considered, making it possible to 
reduce various sources of error to a negligible level. It 
was concluded, in particular, that the most favourable 
orbits are the sun-synchronous, continuous sunlight or- 
bits situated at altitudes between 1390 and 3330 km. 

Since the origination of the SEE concept, the de- 
velopment has focused on critical analyses and critical 
hardware requirements. All indications from this work 
are that the SEE concept is feasible and practicable Q . 
A "blue ribbon" Theory Advisory Group, formed two 
years ago to critique Project-SEE activities and goals, 
has concluded that they are sound. 

This paper presents some new evaluations concern- 
ing the opportunities of the SEE concept and its yet- 
unresolved difficulties. In Sec. 2, for comparison, we 
briefly outline the current status of terrestrial and as- 
tronomical determination of the gravitational interac- 
tion parameters to be measured by the SEE method. In 
Sec. 3, on the basis of computer simulations of Particle 
trajectories, we estimate the requirements for the Shep- 
herd quadrupole moment uncertainty. Sec. 4 shows the 
results of simulations of the measurement procedure it- 
self, which enables us to estimate the possible measure- 
ment accuracy with respect to G and a for A of the 
order of either metres or the Earth's radius. Sec. 5 dis- 
cusses a spurious effect of test body electric charging 
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when the satellite orbit passes through the Van Allen 
radiation belts, rich in high-energy protons. Sec. 6 is a 
conclusion. 

In what follows, the term "orbit" applies to satel- 
lite (or Shepherd) motion around the Earth, while the 
words "trajectory" or "path" apply to Particle motion 
with respect to the Shepherd inside the capsule. 

2. State of the art: a brief survey 

Since gravitational forces are so very small, precision- 
measurement techniques have been at the core of ter- 
restrial gravity research for two centuries. However, 
evidence is increasingly accumulating which indicates 
that terrestrial methods have plateaued in accuracy and 
are unlikely to achieve significant accuracy gains in the 
future J5j. For example, the uncertainty in the gravi- 
tational constant G had been accepted as 128 ppm for 
nearly two decades, and the actual uncertainty in G - 
as indicated by the scatter of results among recent ex- 
periments which claim high accuracy — is roughly the 
same (about 140 ppm). We discuss below the situation 
with respect to several key measurements. 

2.1. Terrestrial determinations of G 

The Luther & Towler determination of G in 1982 ||, 
with the result (6.6726 ± 0.0005) • 10~ n N • m 2 /kg^ 
and other, less precise experiments gave rise to the cur- 
rent official CODATA value of G, viz. 6.67259 • 10~ n 
N • m 2 /kg 2 with an error of 128 ppm. Several still other 
experiments which also claimed high precision were ig- 
nored by CODATA because of inadequate documenta- 
tion of systematic errors. 

There is considerable evidence that the uncertainty 
in G has plateaued at about 100 ppm. At a recent 
(November 23-24, 1998) conference in London, several 
new (1998) determinations of G were reported. The 
obtained values for G (in units of 10 11 N-m 2 /kg 2 ) 
and their estimated error 5G/G in ppm are as follows: 



Fitzgerald and Armstrong 


6.6742 


90 


(New Zealand) § 


6.6746 


134 


Nolting et al. (Zurich) Q 


6.6749 


210 


Meyer et al. (Wuppethal) § 


6.6735 


240 


Karagioz et al. (Moscow) JhJ 


6.6729 


75 


CODATA (1986) 


6.67259 


128 



Obviously, most of the stated errors are of the order 
100 ppm. Moreover, the scatter (1-sigma) about the 
mean is about 140 ppm. Some of the investigators still 
hope for accuracy of 10 ppm. It remains to be seen 
whether they will be able to report error estimates of 
this size or, more importantly, whether their respective 
values for G will actually agree within 10 ppm. 



We note that this analysis is perhaps unduly opti- 
mistic since it excludes one extremely bad outlier: the 
very careful and well documented experiment by the 
Physikalisch-Technische Bundesanstalt in the late 1980s 
and early 1990s, which obtained a series of values for G 
that were consistenly above the results of most experi- 
menters by about 6000 ppm (0.6 % !), while claiming an 
error of about 100 ppm [^3). No explanation for such a 
large discrepancy has been found. 

It might seem that the problems of terrestrial appa- 
ratus must inexorably yield to new technologies — that 
the promise of ever increasing sensitivities would also 
lead to ever improving accuracy. However, this may 
not be true, since it is various systematic errors which 
limit the ultimate attainable accuracy in terrestrial ex- 
periments [^). 

2.2. Terrestrial tests of the equivalence 

principle (EP) and search for Yukawa 
forces 

The EP may be tested by searching for either viola- 
tions of the inverse-square law (ISL) or composition- 
dependent (CD) effects in gravitational free fall. 

In the watershed year of 1986, Fischbach startled 
the physics community by showing that Eotvos's fa- 
mous turn-of-the-century experiment is much less de- 
cisive as a null result than was generally believed Q]. 
Prior to this time, experiments by Dicke jlj| and Bra- 
ginsky [[l6f had demonstrated the universality of free 
fall (UFF) to very high accuracy with respect to several 
metals falling in the gravitational field of the Sun (the 
Eotvos parameter r\ was ultimately found to be smaller 
than 10 -12 ). The interpretation of these results at the 
time was that they validated UFF. 

It was implicit that any violation would have infi- 
nite range, like gravity fl7|| . During the 1970s and early 
1980s there was also a flurry of activity concerning pos- 
sible ISL violations, which eventually led to null results 
at the levels of precision then available (Fujii, fis| , Long 
@) ; 

Since 1986 it has become customary to parametrize 
possible apparent EP violations as if due to a Yukawa 
particle with a Compton wavelength A . This approach 
unites both ISL and CD effects very naturally, while 
the parameter values in the Yukawa potential suggest 
which experimental conditions are required to detect 
the new interaction. 

Following Fischbach's conjecture, ISL and CD tests 
were undertaken by many investigators. Although a 
number of anomalies were initially reported, nearly all 
of these were eventually explained in terms of over- 
looked systematic errors or extreme sensitivity to mod- 
els, while most investigators obtained null results. By 
far the tightest bounds are those obtained by Adel- 
berger and his "Eot-Wash" group at the University of 
Washington 1 20 . This group expects a further improve- 
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ment of at least an order of magnitude |2lJ] . A positive 
result for a deviation from the Newtonian law (ISL) 
was obtained (and interpreted in terms of a Yukawa- 
type potential) in the range of 20 to 500 m by Achilli 
and colleagues p2|; this needs to be verified in other 
independent experiments. 

For reviews of terrestrial searches for non-Newtonian 
gravity sec |l7|, [23], p|. The opportunities of the SEE 
concept in this respect are discussed in Refs. jj], ^| and 
in the present paper. 

The UFF is still in the scope of the current exper- 
imental projects, and the SEE concept suggests here a 
progress of 3 to 4 orders of magnitude as compared with 
Ref. Q. Only one project, STEP (Satellite Test of the 
EP) promises a greater progress but meets some signif- 
icant problems of its own |25| ], connected, in particular, 
with the radiation belts. 

3. Simulations of Particle trajectories 
and the Shepherd quadrupole 
moment 

In the previous studies of the SEE project it was as- 
sumed that the capsule was about 20 m long and the 
initial Shepherd-Particle separation xq along the cap- 
sule axis was as great as 18 m; some estimations were 
also made for 5 m < xq < 10 m. The Shepherd mass 
was taken to be M — 500 kg and the Particle mass 
m = 0.1 kg. The present study retains these values. 

In what follows we describe some characteristic fea- 
tures of Particle trajectories with a goal to determine 
their sensitivity to the uncertainty of the Shepherd 
quadrupole moment J 2 for xq > 5 m. As in our pre- 
vious studies, the capsule diameter is supposed to be 1 
m. 

The reason for considering the quadrupole moment 
uncertainty is technological by origin. Namely it is 
hard to produce a spherically symmetric Shepherd to a 
required accuracy and, instead, it has been suggested 
Q to use a Cook-Marussi stack of cylinders with J 2 > 
, which may be manufactured more easily. A slow 
rotation of the Shepherd with J2 > will stabilize its 
position and orientation. 

The value of J 2 can be provided with some uncer- 
tainty 8J2 ■ To avoid the inclusion of 8J2 in the set of 
parameters to be determined in the experiment, it is 
useful to know which values of 8J2 will be negligible, 
since the growth of the number of parameters leads to 
serious problems in data processing. 

3.1. Equations of motion and the initial data 

Assuming that the relative motion of the test bodies 
inside the capsule occurs in the satellite orbital plane, 
the reduced Lagrangian of the Particle motion reads 



L 



M 



rn 
~2 



f 2 +r 2 (ip + i>f 



G 



G^{l + J 2 (^) 2 p 2 (cose)}(l + 



ae 



-s/A 



(1) 



where (i?, (p) are the Earth-centred polar coordinates 
of the Shepherd in the orbital plane; r — y/ (R+y) 2 +x 2 
and i/j are the Earth-centred polar coordinates of the 
Particle; x and y are the Shepherd-centred Particle co- 
ordinates, where x is the "horizontal" one, i.e., along 
the orbit and simultaneously along the capsule and y 
is the "vertical" one, along the Earth-Shepherd radius 
vector; s — \J x 2 + y 2 is the Particle-Shepherd sepa- 
ration; Mgg , M and m are the Earth, Shepherd and 
Particle masses, respectively; J2 is the quadrupole mo- 
ment of the Shepherd, r s is its radius and P2 is the 
Legendre polynomial 



P 2 (cos( 



3 cos 2 9 - 1 



where 9 is the angle between the line connecting the 
centres of the test bodies and the Shepherd equatorial 
plane. It is easy to see that if the Shepherd sym- 
metry axis is in its orbital plane, then 9 = 9q = 
— arctan(y/x) + <p. If the symmetry axis of Shepherd is 
orthogonal to its orbital plane, then 6 = 0. In general, 
if x is the angle between the Shepherd symmetry axis 
and its orbital plane, then 6 = 9q cos x ■ Hence the in- 
fluence of J 2 on the Particle motion is minimum if the 
Shepherd symmetry axis lies in its orbital plane and is 
maximum if they are mutually orthogonal. 

For simplicity (and taking into account the corre- 
sponding estimate) we neglect the influence of the Par- 
ticle on the Shepherd, so the Shepherd trajectory is 
considered to be given. Then, varying the above La- 
grangian with respect to x and y , taking into account 
that M 3> m and R ^> s, we arrive at the following 
equations of Particle motion with respect to the Shep- 
herd: 



cPx 
~dH 2 



a GM e \ IRCpy 



GM 



x\l + J 2 [ — ) P 2 (cos 



GM . , 
ax — 7T \ 1 + J: 



GM 



x [x(l + 3 cos 29) + 3y sin 20 cos x] 



.,. • ,> 1-2 GM m \ 2Ripx 
— = -2xp + (R + y)W -— jr 



(2) 



d 2 y 
dt 2 



R 



GM 



r. 



yU + J 2 (-) -P 2 (cos 



s J I \ s 

GM (\ 1 
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s' 2 \ s X 



1 + J 2 - P 2 (cos9)\e 
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GMrl t f, ,/a\ 
+ —^-J2 (l + ae s/X J x 

x [3a; sin 26 cos x + V (1 - 3 cos 9)] (3) 

where M = M + m . 

Two kinds of initial conditions for Eqs. (^) and (||) 
were used during the simulations. First, we used the so- 
called "standard" initial conditions, taking the Particle 
velocity components i(0) and y(0) corresponding to 
its unperturbed (i.e., without the M — m interaction) 
orbital motion distinguished from the Shepherd's orbit 
only by its radius (for circular orbits) or major semiaxis 
(for elliptic orbits). Assuming that the Particle motion 
begins right at the moment when the Shepherd passes 
its perigee, these conditions have the form 

x(0) = x 0) y(0) = yo, 

2(1 — ey e (1 — e) 

where oj 2 = GM e /Rq, Rq is the Shepherd orbital ra- 
dius (at the perigee), e is the orbital eccentricity and 

For clearness, the relations (^) are written in the 
linear approximation in the variables x and y . Higher- 
order approximations were used in the simulation pro- 
cess as well. 

The second kind of initial conditions correspond to 
small variations of initial velocities with respect to their 
"standard" values. 

The set of equations (||)-(||) was solved numerically 
using the software developed previously |D to analyze 
the SEE project. 

On the basis of numerical solution of Eqs. and 
, we considered two types of Particle trajectories, 
corresponding to different choices of the initial data: (i) 
approximately U-shaped ones and (ii) cycloidal ones, 
containing loops (see more details on the trajectories 
in Jl|,|l), for orbital altitudes H olb = 500, 1500 and 
3000 km. The uncertainty 8J2 ranged in the interval 
10~ 3 -j- 10~ 5 ; and initial Particle position changed in 
the range 6 m < xq < 18 m, —25 cm < yo < —5 cm. 

For H OI b — 500 km, all U-shaped paths contained 
a sinusoidal component (with the orbital frequency), 
starting at xq < 8 m, while for -ff or b = 1500 and 3000 
km it was present in paths starting at xq < 10 m. 
In other families of U-shaped trajectories a sinusoidal 
component was present only in the case \yo\ > 20 cm. 

3.2. Restrictions on the Shepherd quadrupole 
moment uncertainty 

Small Shepherd quadrupole moment uncertainties SJ2 
create small displacements Sr of a Particle trajectory 
with respect to unperturbed one, r*o(t): 

5r = r. J {t)-f Q {t) 

where fj is the perturbed path. Instead of the full 
displacement 5f, a displacement 5x along the x axis 



Table 1. Displacements of U-shaped trajectories un- 
der 8J2 = 10~ 4 for the Shepherd symmetry axis in its 
orbital plane. The second line shows xq . 



yo 

(cm) 




Sx max x 10 7 


(m) 




18 m 


10 m 


8 m 


6 m 


4 m 


-25 


13.3 


7.52 


7.98 


8.57 


14.90 


-20 


5.51 


4.65 


4.83 


5.51 


5.84 


-15 


1.98 


2.87 


3.03 


3.3 


2.18 


-10 


0.56 


1.59 


1.7 


1.88 


0.66 


- 5 


0.11 


0.59 


0.62 


0.69 


0.14 



Table 2. Displacements of U-shaped trajectories under 
SJ2 = 10 for the Shepherd symmetry axis orthogonal 
to its orbital plane. The second line shows xq . 



yo 

(cm) 




Sx max x 10 7 


(m) 




18 m 


10 m 


8 m 


6 m 


4 m 


-25 


54.2 


62.4 


68.9 


81 


57 


-20 


21.9 


25.8 


28.1 


32.2 


23.3 


-15 


7.94 


10.2 


11.5 


14 


8.7 


-10 


2.25 


3.42 


4.15 


5.45 


2.65 


- 5 


0.45 


0.83 


1.05 


1.47 


0.58 



may be considered since, by numerical simulations, dis- 
placements along the y axis are an order of magnitude 
smaller than Sx. 

Numerical simulations show that in the whole range 
of the above initial conditions the displacement Sx is 
(as it should naturally be) a linear function of 8J2 ; for 
SJ2 = 10~ 4 . For the case when the Shepherd symmetry 
axis is located in the orbital plane, the maximum values 
of Sx for U-shaped trajectories are given in Table 1. 

One can conclude that, if the distance measurement 
error is 10 -6 m, for most of the trajectories the uncer- 
tainty SJ2 = 10~ 4 is admissible. 

The increase of Sx for small values of 1 2/0 1 is ex- 
plained by a large displacement of the turning point 
towards the Shepherd. 

When the Shepherd symmetry axis is orthogonal to 
the orbital plane, these estimations change as shown in 
Table 2. 

For cycloidal trajectories these values are approx- 
imately an order of magnitude smaller than those for 
the U-shaped ones. 

It was also found that Sx decreases with increas- 
ing orbital altitude i/ rb ■ Table 3 shows, as an ex- 
ample, the displacements of U-shaped trajectories with 
x = 18 m for H olh = 500, 1500 and 3000 km and the 
Shepherd symmetry axis located in its orbital plane. 
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Table 3. Displacements of U-shaped trajectories with 
xq = 18 m and different orbital altitudes for 8J2 = 
10~ 4 . The second line shows the values of ff rb • 



2/o 
(cm) 


fernsx x 10 7 (m) 


500 km 


1500 km 


3000 km 


-25 


42.6 


13.3 


4.84 


-20 


11.9 


5.51 


2.18 


-15 


4.02 


1.98 


0.86 


-10 


1.03 


0.56 


0.28 


-5 


0.17 


0.11 


0.07 



The above results show that the Shepherd quad- 
rupole moment uncertainty 8, J 2 may be neglected in 
the SEE experiment with circular Shepherd orbits at 
H or h — 1500 or 3000 km and U-shaped Particle tra- 
jectories if 8J2 < 10 -5 and the position measurement 
error 81 is 10 -6 cm, or 8J2 < 10~ 7 for 81 = 10~ 8 
cm. For cycloidal Particle trajectories or elliptic Shep- 
herd orbits these estimates become 8J2 < 10~ 4 and 
8J2 < 10~ 6 , respectively. For low orbits, -fforb = 500 
km, the restrictions on 8J2 become more stringent: 
8J 2 < 10~ 6 for 81 = 1Q- 6 cm and 8J 2 < 10" 8 for 
81 = 10 -8 cm. However, as is evident from the above 
tables, these requirements may be relaxed by an order 
of magnitude if one discards some trajectories. 

The influence of 8 J 2 on the accuracy of G measure- 
ment may be now estimated as follows. Let some value 
of 8J2 produce the trajectory displacement \8r\ < 8lj 
while the variation SGq of G with the same initial con- 
ditions gives the trajectory displacement \5r\ < 81q- 
Then, keeping in mind the linear dependence of trajec- 
tory displacements on 8J2 and 8G , the accuracy of G 
measurement under the uncertainty 8J2 may be esti- 
mated as 

8G SL SGq 
— < — . 

G SIq G 

Using this inequality and the results of trajec- 
tory simulations, we obtain the following estimates for 
U-shaped Particle trajectories in circular orbits with 
H olh = 1500 km: 

Table 4. Estimates of SG/G in ppm for 8J2 = 10~ 4 , 
when the symmetry axis of the Shepherd lies in (x = 0) 
or is ortogonal to (x = 7r /2) its orbital plane. The 
second line shows xq . 



2/o, 
cm 


x = o 


X = ^/2 


18 m 


6 m 


18 m 


6 m 


-25 


0.88 


0.7 


3.57 


6.7 


-20 


0.27 


0.3 


1.05 


1.73 


-15 


0.06 


0.1 


0.24 


0.44 



One can conclude that the uncertainties 8J2 ^ 10~ 5 
do not create substantial G errors for most of the tra- 
jectories. 



4. Simulations of experimental 
procedures 

This section describes the results of computer simu- 
lations of the whole measurement procedures aimed 
at obtaining the sought-after gravitational interaction 
parameters. These simulations assumed the Shepherd 
mass M = 500 kg, a circular orbit with i? or b = 1500 
km under a spherical gravitational potential of the 
Earth, and a Particle mass of 100 g. Where relevant, 
it is assumed that both the Shepherd and the Parti- 
cle are made of tungsten. Their identical compositions 
are assumed for simplicity since this work is performed 
only for estimation purposes. 

4.1. Equations of motion with Yukawa terms 

We will begin with a presentation of the Particle equa- 
tions of motion in the relevant approximation, includ- 
ing the contributions from hypothetical Yukawa forces, 
taking into account the finite size of the Yukawa field 
sources. 

Let the interaction potential for two elementary 
masses mi and rri2 be described by the potential 



dV 



Yu 



G dm,\dm,2 



-r/\ 



(5) 



where r is the masses' separation, a and A are the 
strength parameter and the range of the Yukawa forces. 
Then for two massive bodies with the radii R\ and R 2 
after integration over their volumes we obtain pfl] 



V 
where 

A = 



Yu Gm 1 m 2 (3if32 - r /\ 
= ae ' 



3 r 



Ri , Ri . , Ri 

— cosh — sinn — 

A A A 



(0) 



(7) 



When A <SC 1 , we have /3j w 1 . This may be the case 
when we consider the interaction between the Shepherd 
and the Particle at a distance of the order of a few 
metres. The radii of the Shepherd and the Particle are 
small: J?i « 18 cm for the Shepherd and R2 ~ 1.1 
cm for the Particle. If the range A is of the order of 
the Earth radius, A « R®, we have /3® = 1.10 and 
/?i.2 = 1 where the indices 1 and 2 label the Shepherd 
and the Particle, respectively. 

The equations of motion are obtained under the fol- 
lowing assumptions. There are two Yukawa interac- 
tions with the parameters Ao and «o referring to the 
Earth- Shepherd and Earth-Particle interactions which 
are the same (due to the assumed identical composi- 
tion for the Shepherd and the Particle) , while A and a 
determine the Shepherd-Particle interaction. The equa- 
tions of motion in the frame of reference connected with 
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the Shepherd, with the same notations x, y, s as pre- 
viously, are 



x + 2uj 2 y + G(mi + mj)- 



3o> 



2 X V 



y , 3w 2 



y — 2lox — 3lu y + G(m\ + 1712)—^ + 



3 1 1 V 



x 
~2 



G(mi 



TO 2 J 



:« 1 



-s/A 



w 2 /? a e- roi/Ao y = 



where u) is the orbital frequency: 



GM a 



'01 



1 + o a a ( 1 + — 

-AO 



''01 



3 — roi/Ao 



(8) 



(9) 



We have neglected the terms quadratic in s/tqi times 
a or ao due to their manifestly small contributions. 

If we set ao = in Eqs. (0), we obtain the equations 
used to describe only the Shepherd-Particle Yukawa 
interaction. One can notice that Yukawa terms are 
roughly proportional to the gradients of the correspond- 
ing Newtonian accelerations, namely, Gmi/s 3 for the 
Shepherd-Particle interaction and GM^/r^ ss lo 2 for 
(say) the Earth-Shepherd interaction. In our case these 
quantities are estimated as 



Gm\ 



2.7-10 



-10 



for s = 5 m, 



8.16- 1CT 7 s~ 2 . 



(10) 



Thus, given the same strength parameter, the Earth's 
Yukawa force is three orders of magnitude greater than 
that between the Shepherd and the Particle, therefore 
one might expect some significant progress in an ISL 
test for A of the order of the Earth's radius. 

The effect of the Earth's Yukawa force is propor- 
tional to the displacements of the statellite along the 
direction of the Earth's radius. Therefore the sen- 
sitivity of the SEE method will increase if one uses 
orbits with eccentricities of the order of 0.01, follow- 
ing Nordtvedt's suggestion j27|. (Larger eccentricities 
would too much disturb the qualitative picture of a 
SEE encounter.) Tentative estimates show that in this 
way one can achieve sensitivities to a ~ 10~ 10 , and 
more thourough studies are in progress. 

Eqs. (||) were used to simulate the measurement pro- 
cedures. 

4.2. Simulations of an experiment for 
measuring G 

The constant G is determined from the best fitting con- 
dition between the "theoretical" (r th (^) = ^ th ) and 
"empirical" (r*;) Particle trajectories near the Shep- 
herd. The fitting quality is evaluated by minimizing a 



Figure 1: Errors 5G estimated by the gradient descent 
( -Rgrad ) and consecutive descent ( _R S ) methods 



functional characterizing a "distance" between the tra- 
jectories. We have considered the following functionals 
for such "distances": 



JV 



i=l 
N 
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The theoretical trajectory depends on the gravitational 
constant G, on the initial coordinates Xq, yo and on the 
initial velocities v x o,v y o- To estimate G, one chooses 
the value for which a "distance" functional in the space 
of the five variables (G,xo,yo, v x q, v v q ) reaches its min- 
imum. 

We carried out a computer simulation of the SEE 
experiment and estimated 5G for a given coordinate 
measurement error (a — 1 ■ 10~ 6 m). As "empirical" 
trajectories, we took computed trajectories, with speci- 
fied values of the above five variables, where a Gaussian 
noise was introduced from a random number genera- 
tor. Independent "empirical trajectories were created 
by non-intersecting random number sequences. The 
functional was minimized using the gradient descent 
method and the consecutive descent method. The start- 
ing value of the "vertical" (along the Earth's radius) 
coordinate, yo , was taken to be 0.25 m, while the hori- 
zontal one, xq , varied between 2 and 18 m. Fig 1 shows 
the dependence of the errors SG/G = R gia d, obtained 
by the gradient descent method and 8G/G = R s , ob- 
tained by the consecutive descent method. All the er- 
rors are estimated by confidence intervals corresponding 
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Figure 2: The SEE method sensitivity to Yukawa forces 
between the Shepherd and the Particle 



to a confidence of 0.95. The mean values of these errors 
are as follows: 

figrad = 4.69 • 1CT 8 , R s = 5.24 -lCT 8 . 

Thus the errors estimated by the gradient and consec- 
utive descent methods are close to each other and are 
about an order of magnitude smaller than the error 
from one-trajectory data. It has been discovered that 
the simulation results strongly depend on the random 
number generator, so that ordinary generators are not 
perfect. 

The use of truncated functionals like (2) has shown 
that a functional incorporating the more informative 
"horizontal" coordinate x leads to estimates close to 
those obtained from the total functional, whereas the 
use of y alone substantially decreases the sensitivity. 
Therefore in practice, to determine G, it is sufficient to 
measure only one of the two coordinates, viz. x . 

Since the "empirical" trajectory is built on the basis 
of a computed one, with a known value of the gravita- 
tional constant Gq , it appears possible to estimate a 
possible systematic error inherent in the data process- 
ing method. The latter has turned out to be in most 
cases much smaller than the random error. This result 
shows the correctness of the methods used. 

As is evident from the results, the best accuracy is 
achieved at values of xo ( ~ the capsule size) about 4-5 
metres. 

4.3. Sensitivity to Yukawa forces with A~lm 

In an experiment for finding a Yukawa interaction be- 
tween the Shepherd and the Particle with the potential 
with /3i t 2 = 1, one computes two theoretical trajec- 
tories: one ignoring the Yukawa forces (x°(ti), y°(ti)) 

and another taking them into account ( x a (ti), y a (ti) ) . 



Figure 3: The SEE method sensitivity to Yukawa forces 
with the range parameter Ao of the order of the Earth's 
radius 



These two computed curves are compared with the em- 
pirical trajectory using the functional Sk (k = 0, a) 
according to ( |Tl"| ) which may be considered as a disper- 
sion characterizing a scatter of the "empirical" coordi- 
nates with respect to the fitting trajectory. This is true 
when the theoretical model is adequate to the real situ- 
ation. In the case k = a the functional Sk = s a has a 
X 2 distribution with n% = 2N — 1 degrees of freedom. 
With k = the parameter a is absent, therefore So 
is distributed according to the x 2 1 &W with Ni — 2N 
degrees of freedom. Then their ratio So/S a = F n2>ni 
will be distributed according to the Fischer law with 
ri2 and n\ degrees of freedom. If an experiment shows 
that, on a given significance level q, the relation 

So/S a > Fni,n 2 ,q (15) 

is valid, one should conclude that a Yukawa force has 
been detected. An equality sign shows a minimum de- 
tectable force on the given significance level q. We 
have assumed q = 0.95. The results of a sensitivity 
computation for different values of the space param- 
eter A are presented in Fig. 2. A maximum sensitivity 
of a = 2.1 ■ 10~ 7 has been observed for A = 1.25 m. 
This value is 3 to 4 orders of magnitude better than the 
sensitivity of terrestrial experiments in the same range. 

These results are based on the measurement method 
which was proposed in the original SEE paper |J; as 
already mentioned, a method involving an eccentric or- 
bit J27J, is much more sensitive and, by our tentative 
estimates, can give an error 6a < 10~ 10 . 

4.4. Sensitivity to Yukawa forces with A ~ Rq 

To estimate the parameter ao in Eqs.(|^), computer 
simulations were carried out using the method as de- 
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scribed above for a, based on the Fischer criterion for 
the significance level 0.95. The range parameter Ao 
varied from (l/32)i?® to 32i?e . Two trajectories with 
the initial Shepherd-Particle separations xq of 2 and 5 
m were calculated. In both cases the impact param- 
eter yo was chosen to be 0.25 m. We used Eqs. (||) 
with a — 0, i.e., excluding the non-Newtonian inter- 
action between the Shepherd and the Particle. As is 
evident from Eqs.(||), the Particle trajectory depends 
on the ratio roi/Ao in the product (roi/Ao) e~ rol / A ° . 
This quantity reaches its maximum at Ao = roi /2 . Our 
calculations have confirmed that a maximum sensitiv- 
ity of the SEE method (3.4 • 10~ 8 for x = 5 m) is 
indeed observed at this value of Ao . This is about an 
order of magnitude better than the estimates obtained 
by other methods. Hopefully this estimate may be fur- 
ther improved by about an order of magnitude by opti- 
misation of the orbital parameters. However, there is a 
factor which can, to a certain extent, spoil these results, 
namely, the uncertainty in the parameter to which, in 
this calculation, was assumed to be precisely known. 

The simulation results are shown in Fig. 3 for two 
trajectories with initial Shepherd-Particle separations 
of 2 and 5 metres. 

5. A possible effect of the Earth's 
radiation belt 

Charged particles, penetrating into the SEE capsule 
from space and captured by the test bodies, create elec- 
trostatic forces that could substantially distort the ex- 
perimental results. Among the sources of such particles 
one should mention (i) cosmic-ray showers, (ii) solar 
flares and (iii) the Earth's radiation belts (Van Allen 
belts). The effect of cosmic-ray showers was estimated 
in Ref. Q and shown to be negligible. Solar flares 
are more or less rare events and, although they create 
very significant charged particle fluxes, sometimes even 
exceeding those in the most dense regions of the radia- 
tion belts, one can assume that the SEE measurements 
(except those of G) are stopped for the period of an in- 
tense flare. On the contrary, the effect of the Van Allen 
belts is permanent as long as the satellite orbit passes, 
at least partially, inside them. 

We will show here that the charging is unaccept- 
ably high at otherwise favourable satellite orbits, so 
that some kind of charge removal technique is neces- 
sary, but this problem may be solved rather easily by 
presently available technology. 

The range of the most favourable SEE orbital alti- 
tudes, roughly 1400 to 3300 km coincides with the 
inner region of the so-called inner radiation belt pq|- 
pl[ , situated presumably near the plane of the mag- 
netic equator. This region is characterized by a consid- 
erable flux of high-energy protons and electrons. For 
a SEE satellite at altitudes near 1500 km the duration 



of the charging periods is about 12 minutes. Maximum 
charging rates occur in the central Atlantic. It should 
be noted that the South Atlantic Anomaly (SAA) — a 
region of intense Van Allen activity which results from 
the low altitude of the Earth's magnetic field lines over 
the South Atlantic Ocean — cannot cause additional 
problems for the SEE experiments. The reason is that 
the SAA mostly contains low-energy protons which can- 
not penetrate into the SEE capsule. 

Electrons are known to be stopped by even a thin 
metallic shell, so only protons are able to induce charges 
on the test bodies. Proton-induced charges on the test 
bodies can create considerable forces. The inner ra- 
diation belt contains protons with energies of 20 to 
800 MeV, and their maximum fluxes at an altitude 
of 3000 km over the equator are as great as about 
3 • 10 cm s for energies E > 10 6 eV and about 
2 ■ 10 4 cm^s" 1 for E > 10 7 eV. At 1500 km altitude 
these numbers are a few times smaller; the fluxes grad- 
ually decrease with growing latitude ip and actually 
vanish at <p ~ 40 . 

It is thus necessary to have some estimates taking 
into account that (i) the capsule walls have a consider- 
able thickness and stop the low-energy part of the pro- 
ton flux and (ii) among the protons that penetrate the 
capsule and hit the Particle, the most energetic ones, 
whose path in the Particle material is longer than the 
Particle diameter, fly it through and hit the capsule 
wall again. As for the Shepherd, its size is large enough 
to stop the overwhelming majority of protons which hit 
it. 

In what follows, we will assume a Shepherd radius 
of 20 cm and a Particle radius of 2 cm and estimate 
the captured charges for some satellite orbits in a cap- 
sule whose walls of aluminium are 2, 4, 6 and 8 cm 
thick. The SEE satellite must actually involve several 
coaxial cylinders for thermal-radiation control, and the 
combined thickness of their walls must amount to sev- 
eral cm. We will assume, in addition, that the Particle 
also consists of aluminium and stops all protons whose 
path is shorter than 4 cm (thus a little overestimating 
the charge since most of protons will cover a smaller 
path through the Particle material). A 100 g Particle 
of aluminium will have a radius of « 2.07 cm. 

It is advisable to determine first which charges (and 
fluxes that create them) might be regarded negligible. 

5.1. Admissible charges 

Let us estimate the Coulomb interaction both between 
the Shepherd and the Particle and between each test 
body and its image in the capsule walls. To estimate the 
spurious effects on the Particle trajectory, it is reason- 
able to calculate its possible displacements due to the 
Coulomb forces from the growing captured charges. We 
assume that the test bodies are discharged by ground- 
ing to the capsule before launching the motion. 
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Criterion. We will call the induced charges, or the 
fields they create, admissible if they cause a displace- 
ment of the Particle with respect to the Shepherd 
smaller than a prescribed coordinate measurement er- 
ror SI (we take here SI — 1CP 6 m) for a prescribed 
measurement time (we take t > 10 4 s). 

A charge on the Shepherd can be estimated as 



(16) 



Qm ~ zSm I J(t, x) dt = eSMF(t, x) 



where e is the elementary charge, x is the capsule wall 
thickness in cm; J(t, x) is the integral proton flux in 
cm~ 2 s _1 after passing through the wall, that is, the 
flux of protons with energies E p > E p (x) where E p (x) 
is such an energy that the proton path in aluminium 
equals x cm; Sm ~ 1256 cm 2 is the Shepherd's cross- 
section; F(t,x) is the fluence, i.e., the total number 
of protons of relevant energies that crosses a square 
centimeter of area for a certain period t . 

In a similar way, the charge captured by the Particle 
may be found as 



< 



eS n 



[J(t, x) - J(t, x+A)] dt 



eS n 



[F(t,x)~ F{t, x+4)] (17) 

where S m « 12.56 cm 2 is the Particle cross-section. 
The subtraction in the square brackets takes into ac- 
count the protons which fly through the Particle with- 
out stopping there. The sign < is used since the effec- 
tive Particle cross-section is smaller than its equatorial 
section. 

The Coulomb acceleration a,Q(t) = qMQm/if m) 
(in the Gaussian system of units) depends on the 
Shepherd-Particle separation r and on the form of the 
function J(t) , which in turn depends on the satellite 
orbital motion. 

The charge- induced Particle displacement is approx- 
imately 



Al 



dt 



dtdQ(t) 



(18) 



since the acceleration is almost unidirectional. If, for 
estimation purposes, we suppose that the flux is time- 
independent, J = Jo = const, and take into account 
that in Eq. @ the difference J(t, x) - J(t, x + 4) 

x) (or even smaller; see particular values in the 
next section), then the resulting displacement is about 



Al 



30 



(19) 



The strong time dependence is explained by the rapid 
growth of the Coulomb force with capturing the charge. 
Numerically, with the above values of Sm and S m , 
taking m = 100 g and r = 1 m (the latter leads to an 
overestimated force since the Particle spends most of 
time at greater distances), one gets: 

J V<0.83-10 18 s 2 cm- 4 . (20) 



For t — 10 4 s an admissible flux is only within 9 
cm s . 

Another undesired effect is that the Particle, being 
charged by the belt protons, will interact with the cap- 
sule walls. This is well approximated as an interaction 
with the Particle's mirror image in the wall, while the 
latter may be roughly imagined as a conducting plane. 
Then, assuming that the Particle is at average at about 
25 cm from the capsule wall and using the same kind 
of reasoning as above, we obtain instead of (EG) 



Jni 4 < 2.07 • 10 19 s 2 cnT 



(21) 



and an admissible proton flux within 45 cm 2 s 1 for 
t = 10 4 s. 

Some more estimates are of interest: if a charge can 
be kept smaller than a certain value, then how great 
may it be to create only negligible displacements? Sup- 
pose that there are constant charges on both the Shep- 
herd {q = <7m) and the Particle (q = q m , m = 100 g), 
then they are admissible according to the above crite- 
rion as long as 



10- 6 CGSE g 2 : 

2 



qNiq m < 2 

q f n < i ■ 10~ 6 CGSE, 



f.io- 



-24 C 2 



(22) 
(23) 



These inequalities follow, respectively, from considering 
the Shepherd-Particle interaction and the interaction 
between the Particle (located at 25 cm from the wall) 
and its image. Thus the maximum admissible Particle 
charge is about 7 • 10~ 4 CGSE, w 1.5 



ing this value, it follows from (22) that the maximum 
Shepherd charge is about 3 • 10^CGSE g rj 5.5 • 10 6 e. 
With these charge values the electric potentials on the 
test body surfaces are 



U M ~ 1-5 • 10~ 4 CGSE g /cm = 45 mV; 
U m w 3.5 • 10~ 4 CGSE 9 /cm = 105 mV. 



(24) 



If by any means the requirements ([22|) , (|23|) are sat- 
isfied (e.g., the potentials are kept smaller than the val- 
ues ([24|)), the electrostatic effect on the Particle trajec- 
tory may be neglected. 

The Shepherd's interaction with its image charge 
induced in its nearest bottom of the SEE experimental 
chamber does not lead to appreciable Particle displace- 
ments. A very demanding requirements on the Shep- 
herd's charge emerges, however, if the SEE satellite is 
used for G-dot determination (whose detailed discus- 
sion is postponed to future papers). One obtains then 



U M < 12 mV. 



(25) 



Evidently, in this case the Shepherd-Particle interaction 
per se is not the determining factor with respect to 
charge limits on the test bodies. 
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Table 5. Average proton fluxes in some satellite orbits at minimum solar activity (x is given in cm, E p in MeV; i is 
the orbit inclination and f2 is its ascension angle, i.e. the longitude at which the satellite crosses the equatorial plane 
moving northward. The notation 1.2(3) means 1.2 • 10 3 , etc. In the first column, the letter 'a' labels equatorial 
orbits, 'b' and 'c' mark less and more favourable orbits (thst is, with greater and smaller numbver of protons), 
respectively, for given altitide and inclination. 
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5.2. Evaluation of charges captured in some 
orbits 

Let us now estimate the charges captured by the Shep- 
herd and the Particle on board a satellite in various 
circular orbits for a single revolution around the Earth, 
a period of about two hours. Actual measurement times 
may exceed this period, but not too much. 

Approximate values of time-averaged proton fluxes 
are presented for some circular orbits in Table 5. 

The fluxes in Table 5 have been obtained using the 
computation software worked out at Nuclear Physics 
Institute (NPI) of Moscow State University, called 
SEE2 (Space Environment Effects 2) and SEREIS 
(Space Environment Radiation Effects Information Sys- 
tem) ||]-& This S0ltware made use OI the NASA 
models AP8-max and AP8-min for calculating the pro- 
ton fluxes [p5| ; however, the latter rest on measure- 
ments performed in the solar maximum of 1970 and 
minimum of 1964, while the NPI software uses some 
modern models of the Earth's magnctosphcre, taking 
into account its evolution on the scale of decades. 

The high-energy particle fluxes in the radiation belts 
are strongly time-dependent; they vary between max- 
ima and minima of solar activity, being, at least at low 
altitudes relevant for a SEE mission, greater at solar 
minima [2-5]. Table 5 shows the fluxes at a solar min- 
imum; similar calculations for a solar maximum show 
smaller values by at average 20-25 per cent; the dif- 
ference exceeds this value only for x — (being thus 
greater for low-energy protons than for higher-energy 
ones) . 

The solar activity varies from one maximum or mini- 
mum to another, the Earth's magnetic field is sensitive 



to all these variations and also varies due to certain 
terrestrial phenomena. Another source of uncertainty, 
probably not a very strong one, is that the shielding ef- 
fect is calculated by SEE2 software for a detector placed 
at the centre of a spherical shell of shielding material, 
whereas the capsule is cylindrical and the angular dis- 
tribution of the proton flow is also uncertain. It is thus 
clear that any values like those presented in Table 5 may 
only serve as a guide, giving correct orders of magni- 
tude. 

The above data make it possible to evaluate the 
captured charges. The results for two orbits, namely, 
1500b and 1500c, are presented in Table 6. These are 
the worst and the best variants of orbits at 1500 km 
among those analyzed (see the caption of Table 5) . 

These and other similar data lead to some conclu- 
sions of importance for the SEE experiments. 

First, the models show zero proton fluxes in equa- 
torial orbits of 500 — 800 km altitudes but indicate 
considerable fluxes at the same altitudes due to cross- 
ing the SAA. It turns out, however, that the SAA is 
overwhelmingly a low-energy phenomenon and almost 
does not affect fluxes on the relevant energy scale be- 
ginning with approximately 65 MeV. Even more, there 
is a very small proton flux due to the SAA even at en- 
ergies over 10 MeV, so that behind a layer of 1 mm the 
SAA influence is already negligible. Therefore, behind 
a thicker metal layer there are actually no secondary 
particles due to SAA protons. 

Second, at 1500 km altitude the fluxes substantially 
depend on the orbit orientation but remain on the same 
scale of a few million protons per cm 2 at energies over 
65 MeV. 



SEE Project: Current Status and New Estimates 



77 



Table 6. Average flux, peak flux and captured charges per revolution in some satellite orbits 
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Third, evidently, at 3000 km altitude both the total 
flux (for x = 0) and especially its high-energy part are 
a few times greater than at 1500 km. 

Fourth, and most important: for all orbits in the 
desirable range of altitudes the charges are quite large 
as compared with their admissible values; they remain 
large even behind rather thick walls. It is thus quite 
necessary to have means to detect and remove the 
charges during the measurements. Moreover, as seen 
from the peak values in Table 6 and from time scans of 
Van Allen charging in orbits of interest (also obtained 
using the above-mantioned software), at a charging 
peak when crossing the magnetic equator the time re- 
quired for the charge on the test bodies to reach its 
maximum allowable values, as listed above, is a matter 
of seconds, not minutes. Therefore the charge must be 
detected and removed as it builds up, on a time scale 
of seconds. 

The detection and measurement of the charge on the 
test bodies can probably be achieved relatively easily 
by an array of minute microvoltmeters attached to the 
inner wall of the experimental chamber. 

Several methods for removing positive charge are 
now being evaluated. A simple and promising method 
may be to shoot electron beams directly at test bod- 
ies. The number of electrons needed is on the order 
of 10 8 /sec. Although this approach has the inherent 
drawback that it requires that an active system must 
perform correctly for many years, it is simple in princi- 
ple and will accomplish the goal. 

6. Conclusion 

Space offers the prospect of quantum leaps in the ac- 
curacy of gravitational experiments. Although space 
is a challenging environment for research, the inherent 
quiet of space can be exploited to make very accurate 
determinations of G and other gravitational parame- 
ters, providing that care is taken to understand the 
many physical phenomena in space which have the po- 
tential to vitiate accuracy. A distinctive feature of a 
SEE mission is its capability to perform such determi- 
nations simultaneously on multiple parameters, making 



it one of the most promising proposals. 

To conclude, we enumerate different SEE tests and 
measurements and show their expected accuracy as cur- 
rently estimated: 

Test/measurement Expected accuracy 

EP/ISL at a few metres 2 • 10~ 7 

EP/CD at a few metres < 10~ 7 (a < 10~ 4 ) 

EP/ISL at ~ _R© < 10~ 10 

EP/CD at ~ i? e < 10~ 16 {a < 10~ 13 ) 

G 3.3 -10- 7 

G/G < 10~ 13 in one year 

The last estimate is only tentative; the subject is 
under study. 
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